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We use multi-frequency information in first year WMAP data to search for the Sunyaev-Zeldovich 
(SZ) effect. WMAP has sufficiently broad frequency coverage to constrain SZ without the addition 
of higher frequency data: the SZ power spectrum amplitude is expected to increase 50% from W to Q 
frequency band. This, in combination with the low noise in WMAP, allows us to strongly constrain 
the SZ contribution. We derive an optimal frequency combination of WMAP cross-spectra to extract 
SZ in the presence of noise, CMB, and radio point sources, which are marginalized over. We find 
that the SZ contribution is less than 2% (95% c.l.) at the first acoustic peak in W band. Under 
the assumption that the removed radio point sources are not correlated with SZ this limit implies 
erg < 1.07 at 95% c.l. We investigate the effect on the cosmological parameters of allowing an SZ 
component. We run Monte Carlo Markov Chains with and without an SZ component and find that 
the addition of SZ does not affect any of the cosmological conclusions. We conclude that SZ does 
not contaminate the WMAP CMB or change cosmological parameters, refuting the recent claims 
that they may be corrupted. 

PACS numbers: 98.65., 98.65.Dx, 98.65. Hb, 98.70.Vc, 98.80., 98.80.Es 



I. INTRODUCTION 

Wilkinson Microwave Anisotropy Probe [1] (WMAP) 
observations of the Cosmic Microwave Background 
(CMB) have ushered in a new era of high precision ob- 
servational cosmology. Such a tremendous increase in 
data quality requires a corresponding increase in the care 
that goes into the data analysis and interpretation. One 
of the lingering concerns surrounding the analysis is the 
residual effect from additional sources of anisotropics. A 
very prominent candidate among these is the Sunyaev- 
Zeldovich effect (SZ). Electrons in hot gas scatter pho- 
tons and distort the blackbody spectrum of the CMB. 
Galaxy clusters, where gas is the hottest, contribute the 
bulk of the effect, shock-heating the gas in their potential 
wells. Radio and microwave band observations pointed 
at known clusters routinely yield SZ detections, so this ef- 
fect is now well established observationally. The scatter- 
ing preferentially raises the energy of the CMB photons, 
but the number of scatterings is low, so the process never 
achieves thermal equilibrium. Therefore SZ appears as a 
CMB temperature decrement at low frequencies and as 
an increment at high frequencies. This frequency depen- 
dence is well known: it is constant in the Raylcigh- Jeans 
(RJ) regime of the blackbody spectrum, and is univer- 
sal (independent of gas temperature or density) for non- 
relativistic electrons. 

In the channels of WMAP, SZ is a temperature decre- 
ment. Table I provides WMAP's frequency bands, and 
the SZ frequency dependence in those bands. The K 
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Band 


v (GHz) 


AT sz /(T CMB y) 


c ; /cf J 


K 


23 


-1.97 


0.972 


Ka 


33 


-1.94 


0.945 


Q 


41 


-1.91 


0.915 


V 


61 


-1.81 


0.820 


W 


94 


-1.56 


0.611 



TABLE I: SZ contribution in the bands of WMAP. We note 
the band name, frequency, temperature perturbation rela- 
tive to the comptonization parameter y, which does not de- 
pend on frequency, and power spectrum relative to the power 
in the Rayleigh- Jeans regime. In the WMAP bands, SZ 
is a temperature decrement. In the Rayleigh- Jeans regime, 
AT sz /(T CMB y) = -2. 



and Ka bands are the most heavily polluted by galac- 
tic contamination, so the best prospects for identifying 
SZ in WMAP are in the differencing assemblies of the 
upper three bands: two assemblies in Q, two in V, and 
four in W. While it is usually argued that SZ is indistin- 
guishable from the CMB in the WMAP channels this is 
actually not so: in CMB temperature units the SZ power 
spectrum increases by 50% from W to Q channel and in 
W it is only 61% of the RJ power. We will show that this 
suffices to place strong constraints on SZ using WMAP 
data alone, a consequence of the remarkably low noise in 
WMAP data. 

In the literature, several groups have attempted to 
identify SZ in WMAP data using cross-correlations with 
other tracers of large scale structure (LSS). In [2] they 
build an SZ template from the XBACs catalog of x-ray 
clusters [3], and fit for the amplitude of SZ, arguing for 
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a 2.5cr detection. In [4] they compute the cross-power 
spectra of the 2MASS extended source catalog [5] with 
WMAP's Q, V, and W bands, then fit for the amplitude 
of the SZ signal, arguing for a 3.1-3.7cr detection of SZ, 
depending on their mask. In [6] they compute the cross- 
power spectrum of the ROSAT diffuse x-ray background 
maps [7] and a WMAP (Q-W) difference map, finding no 
detection. In [8] they compute the cross-correlation func- 
tion of SDSS DR1 [9] galaxy survey data with the V band 
and fit for the amplitude of SZ, finding a 2.7er detection. 
In [10] they build templates for SZ from several optical 
and x-ray cluster surveys, then fit for the amplitude of 
these templates using maps from the W band and [11], 
finding no detection for optical clusters and 2-5a detec- 
tions for x-ray clusters. 

Finally, in [12] they compute the mean temperature in 
concentric rings about APM [13] and ACO [14] groups 
and clusters, noting a decrement they attribute to SZ. 
They interpret this decrement to extend to large scales, 
~ 1 degree, although the covariance of their bins is un- 
clear. Redshift z < 0.2 clusters dominate their sample. 
Their most extreme model assumes that extended SZ 
emission is representative of the temperature and spa- 
tial clustering of gas to z = 0.5. In this model the SZ 
power is 30% of the first acoustic peak of the CMB. Thus 
they conclude the integrity of the WMAP cosmological 
parameter fits may be compromised. 

In this paper we seek to place limits on SZ from WMAP 
data alone, thus avoiding any uncertainties in connecting 
the results of WMAP-LSS cross-correlation to those of 
WMAP auto-correlation. The cross-correlation methods 
used before require a model for relating SZ temperature 
perturbations to some proxies for SZ, which are nonlin- 
ear structures such as clusters, and this connection can 
be quite uncertain. Our method is less model dependent. 
The downside is that WMAP-only methods sacrifice sig- 
nal to noise because they do not focus on matter concen- 
trations, where SZ should be strongest, so our method is 
not optimized for obtaining an SZ detection. However, 
our main goal is to investigate the amount to which the 
CMB analysis may be compromised by the residual SZ 
component, for which we need analysis as model inde- 
pendent as possible. 

Our principal method, described in section II, con- 
structs a linear combination of the WMAP Q, V, and W 
band cross-power spectra which maximizes the contribu- 
tion of SZ, while at the same time minimizing the radio 
point source and CMB contribution. We use this linear 
combination to fit for the SZ power spectrum amplitude, 
using a spectrum shape from halo model calculations [15]. 

In section III we supplement this method by fitting for 
cosmological parameters from the WMAP temperature 
power spectrum using the Markov chain Monte Carlo 
method, allowing for an additional SZ component in the 
power spectrum with a free amplitude to be determined 
by the data. The two methods give consistent results, 



with the latter giving somewhat weaker constraints. Con- 
clusions are presented in section IV. 



II. CROSS-SPECTRA COMBINATION FOR SZ 

We begin our discussion with a simple example. In 
the absence of noise and point sources, the exact SZ 
power spectrum could be computed from the difference 
of any two cross-spectra from different bands. Sup- 
pose we took the differencing assemblies from the Q, 
V and W bands, and compute the cross-power spectra. 
The SZ power is given by, for example, the difference 
(Q1V1); - (W1W2); = 0.25 Cf Z ' RJ in this idealize case, 
where the coefficient is computed from table I. Other 
combinations would also yield the SZ spectrum, with a 
different coefficient. The important things to note is that 
the coefficient is not very small and the CMB cancels ex- 
actly in this combination: any CMB cosmic variation in 
the Q and V bands is the same as in the W band, and is 
subtracted out. In the presence of instrument noise, this 
difference does not give the exact SZ power spectrum, but 
a noisy estimate for it. Given a power spectrum shape, 
we may estimate its amplitude by summing together dif- 
ferent I bins, each weighted to give the proper normal- 
ization and to emphasize bins with high signal to noise. 
Our SZ estimator works in this fashion, except that in 
the presence of point sources, QV— WW is a biased es- 
timate. However, by including other cross-spectra, we 
can form a linear combination of the cross-spectra bins 
to yield an unbiased estimate for the SZ amplitude. In 
the following we find such a combination, and apply it to 
WMAP data. 



A. Estimator and data 

In this section, we introduce our estimator and list the 
data we use. Our estimator, derived in appendix A, is a 
generalized version of the point source estimator of [16]. 

We want to generate an estimate for SZ from the 
WMAP cross-spectra. We postulate that the data is 
the sum of the CMB and two contaminants, radio point 
sources and SZ. We marginalize point sources and esti- 
mate SZ in our main case, which we illustrate in detail 
below. We consider other cases also, but for brevity omit 
their details, which are similar. We can write the cross- 
spectra as 

(cf) = cf CMB + c ; Mrc + c;< sz , (i) 

showing explicitly the contribution from each part of the 
signal. Here the multipole bin is denoted by I and the 
pair of differencing assemblies in the cross-correlation is 
denoted by i = — W1W2, Q1V1, etc. No auto-power 
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spectra are included, so we need not worry about noise 
subtraction. 

We will marginalize over the CMB spectrum, which 
we denote by Cf MB . The window functions for each 
differencing assembly pair are w = {w^,}. We boldface 
w because later we will think of it as a matrix. Therefore 
the contribution to the cross-spectrum from the CMB is 



ct CMB = £4<^ MB 



(2) 



The spectra we use are in temperature units, and have 
already been beam-deconvolved, so the window functions 
w ii' = are trivial. However, it is necessary to keep 
the window functions explicit in the manipulations that 
follow. 

We denote the amplitude of the point source power 
spectrum by A. This amplitude relates to the cross- 
spectra via the frequency and shape dependence S = 
{SI}. Later we will think of S as a vector. We take radio 
sources to have a white noise spectrum with power law 
frequency dependence, given by: 



C 



Si A 



S) = ( II! ] (12! 



(3) 



where cross-spectrum i has channels at v il and v i2 . The 
units of A are temperature squared. Well-resolved point 
sources have already been masked from the maps before 
the evaluation of the cross-spectra, so A represents un- 
resolved sources only. [2] found f3 = —2.0 and = 45 
GHz for the resolved sources. Following [16] we take the 
same for the unresolved sources. 

We describe amplitude of SZ with B 7 the ratio of the 
SZ power spectrum and the predicted spectrum, assum- 
ing they have the same shape. Thus the SZ amplitude 
is dimensionless, and has a theoretically predicted value 
B = 1 for <r 8 = 0.9 using the halo models of [15]. In 
our notation, this amplitude also relates to the power 
spectra by the frequency dependence and shape for SZ, 
which we label Z = {Z\}. The frequency dependence of 
a temperature perturbation due to SZ is 



Ar sz 

^cmb « - 2 /( x )' 



(4) 



where the frequency dependence relative to the RJ regime 
is given by 

x/2 



f{x) = 2 - 



(5) 



tanh(a;/2) 

where x = hu/ksTcMB- Note / — > 1 in the RJ limit 
x — > 0. Thus the SZ contribution to the cross-spectrum 
is: 



£,i,SZ 



Z\B 



where the spectrum shape C ; SZ ' RJ i s from the halo model 
prediction in the RJ regime using <r 8 = 0.9, as shown 
in Figure 3. The shape of SZ is roughly (7 ; SZ - RJ cx l^ 1 , 
although it becomes slightly steeper for I greater than a 
few hundred. 

We organize the binned cross-spectra C\ into a data 
vector D = {C/}. We use a Gaussian model for the 
likelihood C of the power spectrum: 



-21og£ oc [D — (D>] f IT 1 [D - (D) 



(7) 



Z) = C\ 



SZ.RJ 
I 



/0ii)/(^ 2 ). 



(6) 



where the covariance £ = ((D — (D))(D — (D))t) can 
be written as £ = {£",}. We derive B, an unbiased 
estimator for B, and its covariance S B in appendix A. 
The main result is: 

B = (Z t FZ)" 1 Z t FD 
S B = (Z+FZ)- 1 , (8) 

where we have defined the auxiliary matrices 

F = E-ES(S t ES)~ 1 S t E 

E = XT 1 - (w t E- 1 w)" 1 w t ir 1 (9) 

In this notation, we consider D, S, and Z as column 
vectors with single index il, and w as matrix with in- 
dices il and S, E, and F are matrices with indices 
M and \'V . This estimator marginalizes CMB and point 
sources, which is the most conservative treatment, since 
it assumes we know nothing about these components. 
A more aggressive treatment would be to estimate all 3 
components simultaneously, but we defer this approach 
to a future analysis. Note that B is a linear combination 
of the cross-spectra D, with weights (Z^FZ^Z^F. 

The data we use consist of the 28 cross-power spec- 
tra from the eight differencing assemblies in the WMAP 
Q, V, and W bands. These spectra are provided at the 
Legacy Archive for Microwave Background Data Analy- 
sis [17]. A galactic foreground model [2] has already been 
subtracted out. WMAP's temperature power spectrum 
is a linear combination of these 28 cross-spectra [16]. 

We bin the spectra in /, accounting for the number of 
modes at each I. The Legacy Archive does not provide the 
cross-spectrum covariance, which we need for £ in equa- 
tion (9), so we estimate the covariance from the data. We 
assume the covariance is diagonal in I, which is a good 
approximation [16]. Then in a single bin, we use the 
dispersion of the cross-spectra about WMAP's combined 
temperature power spectrum to estimate £ for that bin. 
(For this purpose we must first un-correct the combined 
spectrum for the point source contribution.) This pro- 
cedure to obtain the covariance works best if the power 
spectrum variance does not change much within a bin 
and if the bin contains enough Vs to get a low- noise esti- 
mate of the variance. Our bin width of Al = 40 is fairly 
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FIG. 1: Reconstruction of the SZ power spectrum in the RJ 
regime, in bins of Al = 40 (red error bars). There is no 
visible SZ detection, which would have been positive for this 
combination. For comparison, we have included the WMAP 
combined temperature power spectrum, binned the same way, 
the WMAP best-fit theoretical CMB power spectrum, and a 
theoretical prediction of SZ in the RJ regime for as = 0.9 
(blue dotted), calculated from a halo model. 



narrow, and gives a covariance which is numerically sta- 
ble in our subsequent calculations. This technique does 
not account for the cosmic variance of the CMB, but as 
we note in the derivation in the appendix, our estimator 
is insensitive to this source of variance. This makes sense 
because the CMB is completely projected out, as in our 
QV— WW example at the beginning of section II. 

As a test of our estimator, we repeat the estimate 
of [16] of the power spectrum amplitude of unresolved 
point sources. In this case we neither marginalize nor 
estimate SZ, but assume B = in equation (6). We 
find a point source amplitude of A = 0.016 ± 0.001 /xK 2 , 
which is roughly consistent with their quoted value of 
A = 0.015 ^K 2 . This gives us confidence in our estima- 
tor, despite our covariance matrix constructed from the 
data. 



B. Results 

We show our estimates for SZ on a bin-by-bin basis, be- 
fore combining different multipole bins to improve statis- 
tics. We marginalize over the CMB and the point source 
amplitude. Our best estimate for the binned SZ power 
spectrum is shown in Figure 1. Immediately we can see 
that the data do not tolerate a large SZ contribution. 

Next we turn to our main case, where we combine all 
the bins together and estimate the SZ amplitude. Fig- 
ure 2 shows the weight (Z' ! FZ)~ 1 Zi' f F for each cross- 



spectrum in the amplitude estimator. The weights have 
been divided into groups based on the bands. For plot- 
ting, we sum the weights in each group. The bulk of 
the weight comes from I = 100 to I = 400. At lower I 
there are fewer modes and the SZ contribution is low. At 
higher I the statistics are limited by detector noise. 

The weights are difficult to interpret heuristically. 
From the frequency dependence, the strength of SZ de- 
creases in order of QQ, QV, VV, QW, VW, and WW. 
We would expect our combination of cross-spectra would 
have the SZ-stronger bands minus the SZ-weaker bands. 
So it is easy to understand the positive weights of QV and 
QW and the negative weights of VW and WW. Point 
sources are strongest in Q-band, so a negative weight 
for QQ makes sense in terms of a point source correc- 
tion. QW and VW are positive except for bins which are 
strong in QV. These negative bins may also represent a 
point source correction 

For the amplitude of SZ relative to the expected RJ 
amplitude from the halo model with WMAP parameters, 
our estimator gives us B = —0.042 ± 1.685. The error 
is large, and is consistent with both no SZ and the ex- 
pected B = 1. Including only the physical region B > 0, 
we integrate the likelihood until we include 95% of the 
probability. From this we quote an upper limit of 

B < 3.3 (95%). (10) 

We plot this limit, along with some models of the SZ 
power spectrum in Figure 3. At the first acoustic peak, 
SZ is less than 3% of the CMB in the RJ regime and less 
than 2% of CMB in W band, on which cosmological pa- 
rameter estimation is heavily based. We find no evidence 
for a large SZ contribution. Assuming Cf z oc (as) 7 , this 
gives a limit of <j$ < 1-07 at 95% confidence. To this one 
should add an additional modeling uncertainty at the 
level of 10% based on the comparison of predictions from 
different simulations [15]. There is a caveat in the up- 
per limit derived here in that we are working with power 
spectra based on masked maps in WMAP, with more 
than 200 radio point sources removed. If the SZ signal 
is correlated with these point sources, which could hap- 
pen if these radio sources sit in massive clusters [18], then 
more SZ may have been removed than expected based on 
the sky fraction of the mask. This would only affect the 
us limits and not the SZ contamination on the primary 
CMB in the WMAP power spectra. Since we are mostly 
concerned with the latter we do not explore this issue 
further. Note that the upper limit on a$ is already com- 
parable to the predicted value based on detections from 
CBI and BIMA, which gives ct 8 - 0.95 - 1.05 [15, 19]. 
It is remarkable that WMAP first year data have suffi- 
cient sensitivity to place constraints on the SZ amplitude 
comparable to other small scale surveys. 

In the next application we jointly estimate the point 
source and SZ amplitude, rather than directly marginal- 
izing over the point sources. We find the two parameters 
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FIG. 2: The contribution to the SZ amplitude estimator B 
from the cross-spectra. The horizontal axis is the multipole I. 
The vertical axis show the weight of each spectrum in fiK~ 2 . 
We have grouped the cross-spectra by band as noted. Within 
each group, we have summed the weights of the individual 
spectra. In the top frame, we have included the CMB power 
spectrum (arbitrary units) to give scale for I. 




FIG. 3: Limit on SZ power spectrum. We present a 95% 
upper limit on the amplitude of the SZ power spectrum from 
a combination of WMAP band cross-spectra. We show the 
temperature fluctuation in the Rayleigh- Jeans regime. For 
comparison we display the WMAP best fit power spectrum 
and SZ power spectra from a halo model calculation using 
cr 8 = 0.9, 1, 1.1. 

to be somewhat correlated (Figure 4), but not to the 
point of allowing very large SZ amplitude. An extremely 
strong SZ power is not allowed. 



III. MARKOV CHAIN ANALYSIS 

In this section we estimate how much SZ signal is in the 
WMAP Ci power spectrum, and investigate the effect of 
SZ on the determination of the cosmological parameters. 
For this purpose we use the Markov Chain Monte Carlo 
(MCMC) approach, using software described in more de- 
tail elsewhere [20, 21]. 

We ran two MCMCs, one without SZ and one allow- 
ing for an unconstrained SZ contribution. We built a 
third chain from the second by importance sampling, al- 
lowing for an SZ component but constraining it to lim- 
its derived based on frequency information in the pre- 
vious section. We used the WMAP likelihood routine 
[16, 22]. Each of the chains contains 100,000 total chain 
elements. The success rate is 45-55 percent, the corre- 
lation length is 13-20 elements, and the effective length 
is 5,000-10,000 elements. Each chain comprises 23 inde- 
pendent sub-chains and, in terms of Gelman and Rubin 
i?-statistics [23], we find the chains are sufficiently con- 
verged and mixed (R < 1.01, compared to the recom- 
mended value of R < 1.2). 

In the second chain, we added to the power spectrum 
an SZ-shaped contribution, parameterized in terms of 
amplitude B (equation 6). The WMAP power spectrum 
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FIG. 4: A joint estimation of the SZ power spectrum ampli- 
tude and the point source spectrum amplitude. The ellipse 
encloses 68% of the likelihood. The best fit point is shown 
by a "+". The vertical line at B=l shows the theoretically 
predicted SZ amplitude. The horizontal line shows WMAP's 
value for point sources, which does not take into account SZ. 
Because the power spectrum is a positive quantity, note that 
portions of the plane where either axis is negative are unphys- 
ical. 

combines the Q, V and W bands in different ratios at 
each I, so the shape of the SZ contribution to the WMAP 
power spectrum is not exactly given by the single fre- 
quency SZ template, because the effective frequency for 
every I varies. This dependence is small and we ignore 
it here. We find that the SZ contribution to the WMAP 
combined temperature power spectrum, dominated by V 
and W channels, may be approximated as 75% the contri- 
bution in RJ. One could also add additional CMB exper- 
iments (e.g. CBI, VSA, etc.) into the analysis, but this 
would incur complications to account for the different 
frequencies of these experiments. In the third chain we 
add our multi-frequency analysis limit as an additional 
constraint. 

We consider only the simplest model required by the 
data plus the SZ component, since we want to analyze 
the effect of the latter on the cosmological parameters. 
WMAP temperature data require neither tensor modes 
nor curvature nor running of the primordial power spec- 
trum of the scalar perturbations, so we do not consider 
them. 

We work in a seven parameter space: 

p = {uj b , w cdm , fi m , r, A s , n s , B}. (11) 

Here u>b — Qbh 2 is the baryonic content of the universe, 
w c dm = n c d m /i 2 is the physical density of the cold dark 
matter content, f2 m = f2 c dm + f2& — 1 — is the matter 
density today, r is the optical depth to reionization, A s is 
the amplitude of the primordial scalar perturbations, n a 




0.02 0.022 0.024 0.026 0.028 0.03 0.032 



FIG. 5: We show two-dimensional contours of 68% and 95% 
confidence levels in (u>b,B) plane. The larger (blue) contours 
show the degeneracy using the WMAP combined temperature 
power spectrum without any prior on SZ. The smaller (red) 
contours include our multi-frequency cross-spectrum analysis 
as a prior. B = 1 is expected from halo models for SZ when 
cr 8 = 0.9. 



is the primordial slope, and as before B is the amplitude 
of the SZ power spectrum. 

To reduce the degeneracies while running the MCMCs, 
we use ojb, w c dm, angular size of the sound horizon O s , 
lgA s , n s , lgA s - t - 0.51g(w fc + Wcdm), and B, instead 
of the parameters in equation 11. We adopt broad flat 
priors on these parameters, and additionally require r < 
0.3. 

We find that the amplitude of an SZ-shaped compo- 
nent to the WMAP power spectrum is limited to B < 7.1 
at 95 percent confidence (see Figure 3). This limit means 
that the contribution to the WMAP temperature power 
spectrum at the first peak is below 5 percent. This is a 
weaker limit than from the combination of cross-spectra, 
where a factor of 2 better limit was found (equation 10). 

Table II shows the comparison of the two MCMCs, 
showing the effect of including SZ in the analysis. The in- 
clusion of the additional parameter B into the likelihood 
analysis affects only the determination of the baryon 
physical density Uf,. Without SZ we find ujb — 0.0235, 
whereas with SZ we find Ub = 0.0251, which is shifted 
by about 1.5cr away from the earlier value. The confi- 
dence contours in the (u>b, B) plane are shown in Fig- 
ure 5, showing that there is a degeneracy between these 
two parameters. However, we can and should also use the 
constraint from our multi-frequency cross-spectrum anal- 
ysis as a prior. We can include the Gaussian likelihood 
for B from section II, and perform importance sampling 
of the chain with SZ. We then find uj b = 0.0243, different 
from the case without SZ by 0.6a. The likelihood con- 
tours including the prior are also shown in Figure 5. The 
other parameters are much less affected by the SZ. The 
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no SZ with SZ, no SZ prior with SZ prior 

B < 7.1 (95%) < 2.9 (95%) 

, , v 1n 2 9 or+0.14 +0.28 9 ri +0.21 +0.47 9 4 o+0.15 +0.30 

w t> A ±u ^- OJ -0.13 -0.26 ■ i - JX -0.18 -0.33 ■ i -^°-0.15 -0.29 

O n 9/l r+0.07 +0.15 n 9O/.+0.07 +0.15 n 9/1 O+0.07 +0.15 

iL m u - zq:d -0.06 -0.10 u - za ^-0.06 -0.10 u - z ^°-0.06 -0.10 

n i i i +0.016 +0.033 n i -i i +0.016 +0.033 n i -i i +0.016 +0.033 

w cdm U - 1 - L1 -0.015 -0.029 U -- L11 -0.016 -0.030 U -- L11 -0.016 -0.031 

_ n 1Q+0.07 +0.10 n 9n +0.07 +0.10 n i q+0.07 +0.10 

' u ' ±a -0.08 -0.14 u - zu -0.08 -0.14 u - ly -0.08 -0.14 

_ n ss+ ' 12 + - 25 n «R+ - 12 + - 24 n 07+0.12 +0.23 

a s U.88_ 011 _ 20 u.8b_ 011 _ 22 \J.Hi_ 011 _ 22 

1. n 7/1+0-06 +°- 12 n 7C+0-07 +0.14 n 7c -+0.06 +0.13 

" u - ,y: -0.05 -0.09 u -' u -0.06 -0.11 u -' d -0.05 -0.10 

r> QQ+0.04 +0.07 n qq+0.04 +0.07 n qq+0.04 +0.07 

"s u - 3a -0.04 -0.07 u - aa -0.04 -0.07 u - aa -0.04 -0.07 



TABLE II: The first two columns contain the median value and 1- and 2a constraints on cosmological parameters for two 
MCMCs without and with a Sunyaev-Zeldovich component in the C; power spectrum from WMAP data alone. For both chains 
there was an imposed prior of r < 0.3. The third column shows the constraints when the limit on SZ from our cross-spectrum 
estimator is applied as a prior. 



next largest deviation from the first chain to the second 
is a shift in the median h of about 0.3a towards higher 
value, which is not statistically significant and the effect 
is even smaller if the constraints from previous section 
are added to the chains. 



IV. CONCLUSIONS 

We have analyzed the WMAP power spectra in two 
ways and found them to be free of serious contamination 
from SZ. The more powerful and less model dependent of 
the two methods is to use multi-frequency information: 
by combining the WMAP cross-spectra, we limit at 95% 
confidence the amplitude of SZ to be below 2% of the 
CMB at the position of first peak in W band. By search- 
ing for an SZ-shaped component in the WMAP combined 
spectrum with a Markov chain Monte Carlo we also do 
not find any evidence of a signal, but we can only set 
a weaker limit. Combining the analyses, we show that 
the cosmological parameters are not affected by the SZ 
within the range allowed by the multi-frequency analysis. 

There are several possible improvements in the analysis 
that could further tighten the limits and may be applied 
to the second year data once they become available. In 
addition to having lower noise, the two year data are 
also expected to have better control of systematics such 
as beam uncertainties and noise correlations. 

First, we would like to improve our method for obtain- 
ing the cross-spectrum covariance matrix. The estimate 
for the covariance is noisy, which makes our weights for 
combining the power spectra noisier than they should be. 
If we double the size of the bins we average over to ob- 
tain the covariance, our limit changes slightly (B < 4.3 
opposed to B < 3.3 at 95%) and the point source es- 
timate is virtually unchanged. However, it is impossi- 
ble to tell if this effect is due to a covariance matrix 



with less noise (which is better) or due to weights whose 
wide bins combine high- and low-signal-to-noise multi- 
poles less optimally (which is worse). This can be im- 
proved if the WMAP second year data release includes 
the cross-spectrum covariance matrix. 

The second possible improvement of our analysis is to 
improve the frequency dependence of point sources or the 
WMAP beam deconvolution, both of which can bias our 
estimate. Our estimator is only unbiased if the models 
one applies to it are faithful to the data. The models we 
have used are the best available, and this information will 
improve in the second year data. It seems unlikely that 
these improvements will strongly modify the conclusions 
in this paper. 

The third improvement is to perform the analysis on 
the maps which do not have point sources excised. As 
discussed above, masking of resolved point sources from 
the WMAP data may reduce the SZ signal if the two are 
correlated. This would weaken our limit on cr 8 from the 
absence of SZ signal, but does not change our conclusions 
about the effect of SZ on the CMB and cosmological pa- 
rameters. Still, given that the current limits are reaching 
the levels of interest it is worth exploring this possibility 
further. 

The fourth improvement is to verify the assumed shape 
for the SZ power spectrum. This is currently based on 
halo models, since hydrodynamic simulations cannot sim- 
ulate sufficiently large volumes to make predictions re- 
liable. If the power spectrum has a radically different 
shape, our limit is hard to interpret, but the level of 
contamination on primary CMB is likely to remain un- 
changed, as is clear from Figure 1. As an example we 
estimated the SZ in bins of Al — 120. In the bin con- 
taining the first acoustic peak, we limit the SZ power 
to 4.9 times the halo model power predictions at 95% 
confidence or 4% contamination in W. It is clear that as 
long as the SZ spectrum is smooth the limits on the con- 
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tamination are not going to change significantly over the 
range with best sensitivity at 100 < I < 300. 

Finally, our procedure to marginalize over CMB and 
point sources is conservative and the errors could be fur- 
ther reduced if a joint estimation is used instead. We 
know something more about the CMB power spectrum 
than just its frequency dependence, and in this case we 
can use this information. We also note that a relatively 
trivial modification to WMAP's method for estimating 
the power spectrum could eliminate any contaminant 
with an SZ-like shape and frequency dependence, much 
in the same way that they eliminate point sources. How- 
ever, the fact that we find no residual SZ effects implies 
that this procedure is not really required and justifies the 
treatment of the WMAP team in ignoring SZ. 

Even with future data detecting the SZ power spec- 
trum from WMAP alone may be difficult. If the bulk 
of the covariance is from noise, it would take WMAP 
roughly 3 years to get a la detection of SZ using the 
cross-spectrum combination method of this paper, if 
as = 0.9. Elimination of non-noise sources of covariance, 
such as uncertainties remaining in the beams, may help. 
If us = 1 then a 2a detection should be possible with 4 



year data and one should start seeing a hint of the signal 
already with 2 year data. Theoretical predictions remain 
somewhat uncertain, so it is possible that these predic- 
tions have a systematic uncertainty of around 10%, but 
it is clear that it will be worth looking for SZ in future 
WMAP releases. 
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APPENDIX A: CROSS-SPECTRUM 
CONTAMINANT ESTIMATOR 

In this appendix we present a generalized version of the 
point source estimator of [16], using a vector notation. 
We discuss contaminants in general and do not mention 
point sources or SZ specifically. We assume we know the 
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power spectrum shape and frequency dependence of all 
contaminants. 

Our data are the cross-power spectra from the experi- 
ment. Let vector D = {C\} be these spectra. The mul- 
tipole (or multipole bin) is denoted by I and the cross 
correlation is denoted by i = W1W2, W1V1, etc. We 
use a Gaussian model for the likelihood £(D) of the set 
of cross-spectra: 

-21og£oc[D-(D)] t S- 1 [D-(D)]. (Al) 

where the covariance £ = ((D — (D))(D — (D}) t ) can be 
written asS= {£",}. 

We postulate that the data D is the sum of the CMB 
and a number of contaminants. For each of these contri- 
butions to the data we need a model, and parameters for 
that model. Some of these parameters we will estimate, 
and some we will marginalize. 

First we consider the CMB. The parameters which de- 
scribe the CMB are the band powers: Ci. We organize 
these band powers into a vector C = {C;}. We will 
later marginalize over these parameters. To relate the 
CMB parameters to the data we use the window func- 
tion w = {w^,} for each cross spectrum channel pair. 
We may think of the window function as the response of 
the instrument to a given set of CMB C/'s. In the ab- 
sence of noise and contaminants, the data would be given 
by the matrix multiplication D = wC. 

We use as similar notation for the contaminants. Wc 
describe each contaminant by a power spectrum shape, 
a frequency dependence, and an amplitude. We take the 
shapes and frequency dependences as given, and the am- 
plitudes as parameters. We divide the amplitude param- 
eters into two groups, those to marginalize and those to 
estimate. The amplitude parameters to marginalize wc 
denote by vector A = where a runs over the com- 

ponents to marginalize. The amplitude parameters to 
estimate we denote by vector B = {Bp}, where (3 runs 
over the components to estimate. 

We may think of the shape and frequency dependences 
as the response of the instrument to a given set of con- 
taminant amplitudes. Thus the shape dependence al- 
ready includes the influence of the window function. Sim- 
ilarly to the amplitudes, we divide the shape and fre- 
quency dependences into two groups. The shape and 
frequency dependence of the contaminants to marginal- 
ize we write as S = {Sl a }. The shape and frequency 
dependence of the contaminants to estimate we write as 
Z = {Z\p\- In the absence of CMB and noise, the data 
would be given by D = SA + ZB. 

If we include all components and uncorrelated noise, 
then we can write the expected data as 

(D) = wC + SA + ZB. (A2) 

We are considering cross-spectra only, and have included 
no noise term. 



We substitute into the likelihood expression: 

-21og£ oc [D - (wC + SA + ZB)]* 

[D - (wC + SA + ZB)] . (A3) 

We seek to estimate B. This we accomplish by repeat- 
edly completing the square in the likelihood expression, 
as follows. Note that completing the square to marginal- 
ize out a component is equivalent to letting the covari- 
ance for that component tend to infinity [24]. Thus we 
disregard all information about that component. The 
likelihood may be rewritten as a quadratic equation in 
the CMB power spectrum, 

-21og£ oc CtaC + l(b + C + C+b) + c^^c, (A4) 

where we have introduced the shorthand 
a = w ji 'E~ 1 w 

b = -2w t S- 1 (D - SA - ZB) 

c = D-SA-ZB. (A5) 

We complete the square for variable C: 







t 




2 log£ oc 


C + ia-b 


a 


C + ia-b 











+c t 5T 1 c--b t a -1 b (A6) 
4 

Note that the final term may be rewritten in terms of our 
auxiliary variable c, 

Jbta^b = c+ (S- 1 ™ (w^S^w) _1 ti^ST 1 ) c. (A7) 

Thus if we define the symmetric matrix 

E^S" 1 -■£- 1 w(w^- 1 w)~ 1 w'f£- 1 , (A8) 

we may compactly express the log likelihood as a term 
which depends on the CMB C, and a term which docs 
not. 



C + ia-b 


t 

a 


CH 


■H 











(A9) 

Let us define a marginalized likelihood, Cc = J dC C. 
Integrating over all C, we find 

-21og£ c oc c f Ec 

oc [D-SA-ZB]tE 

[D-SA-ZB]. (A10) 

Note the similarity to equation (A3). The matrix E is 
the new inverse covariance, once the CMB is marginalized 
out. 
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We wish to repeat this sequence to marginalize the 
variable A. Thus we write 

-21og£ c oc AtfA+ i(gtA + Atg) +h t Eh (All) 

where we have introduced 
f = S^ES 

g = -2S t E(D-ZB) 

h = D-ZB. (A12) 

Again we complete the square, now for variable A: 

-21og£ c oc [A + f-^f [A + f-ig] 

+htEh- igtf-ig. (A13) 

Noting the last term may be re-written in terms of h, 

Jg f f _1 g = h f (ES (S^S)" 1 S f E) h, (A14) 

we write with analogy to equation (A8), 

F = E ES (S f ES) _1 S f E. (A15) 

Define Cc,a = J dCdA C = J dA Cc, and we have 
integrated out our contaminants: 

-21og£ c ,A oc h+Fh 

oc [D - ZB] f F [D - ZB]. (A16) 

To estimate the amplitudes B, we express the 
marginalized likelihood as 

-21og£ c ,A oc BtpB+ i(qtB + Btq)+r (A17) 

where we have introduced 

p = Z f FZ 
q = 2Z f FD 

r = D^FD. (A18) 
We complete the square one final time for variable B: 

-21og£ c ,A cx [B-ip-iqjtptB-^p-'qRr-^qV'q- 

(A19) 

Now£ c ,a(B) occxp[-i(B-(B))t(S B )- 1 (B-(B))]. If 
we estimate 

B = (Z t FZ)" 1 Z t FD, (A20) 



we have (B) = (B), which indicates our estimator is 
unbiased. Moreover, the estimator is unbiased even if 
the original covariancc S from equation (Al) is incorrect. 
This is shown by integrating the (flawed) estimator with 
the likelihood to take the ensemble average. We also note 
the covariance on the B parameter estimates, 

S B = (ZtFZ)- 1 . (A21) 



We note that neither the estimator B nor its variance 
S B are sensitive to the cosmic variance from the CMB. 
This makes sense because the CMB is completely pro- 
jected out. The estimator does not care about the value 
of any CMB multipole, so the cosmic variance of the mul- 
tipoles is also immaterial. 

This independence from cosmic variance may be shown 
algebraically. The estimate appears to depend on the cos- 
mic variance contribution to E (equation A8) through F 
(equation A15). However, if we explicitly write the cos- 
mic variance S c = <(C - (C))(C - (C))t) contribution 
to the variance, we see that it projects out of the esti- 
mate. Let us define £' as the balance of the variance, 
the part not due to cosmic variance of the CMB. Write 
the covariance: 

£ = w^ c w^ + (A22) 
then substitute into E, 

e = (tus^+s'r 1 

™t (w^w^ + . (A23) 

If we were to expand each of the matrix inverses in geo- 
metric series, we would find that E does not depend on 
S c at all: 

e = (E'r 1 - (st 1 w (™ f (s'r 1 ™) _1 ^ (st 1 . 

(A24) 

We see that the cosmic variance of the CMB has been 
projected out of the estimator. Thus it cannot impact 
the estimate B, or its variance S B . 



